AD- A 144  336  PATH  INTEGRATED  GROWTH  OF  AURORAL  KHOMETRIC  RADIATION  1/| 
IU)  IOWA  UN IV  IOWA  CITYVDEPT  OF  PHYSICS  AND  ASTRONOMY 
N  OMID1  ET  AL.  FEB  84  U.  OF  IOWA-84-9  NOOO 14-76-C-00 16 

F/G  20/14 


UNCLASSIFIED 


NL 


AD- A 144  336 


PATH  INTEGRATED  GROWTH  OF 


AURORAL  KILOMETRIC  RADIATION 


N.  Omidi  and  D.  A.  Gurnett 


UVVERSITY  cu 

IWW 

A— <  7X 
°0 (SIDED  t6h 


U.  tftlowa  84.9 


Department  of  Physics  and  Astronomy 

THE  UNIVERSITY  OF  IOWA 


Iowa  City,  Iowa  52242 


Vhi;  document  has  been  apt 
fLi  public  release  and  sale;  i 
distribution  rs  unlimited. 


84  08  09  052 


U.  of  Iowa  84-9 


PATH  INTEGRATED  GROWTH  OF 
AURORAL  KILOMETRIC  RADIATION 
by 

N.  Omidl  and  D.  A.  Gurnett 


Department  of  Physics  and  Astronomy 
The  University  of  Iowa 
Iowa  City,  IA  52242 


Submitted  for  publication  to  Journal  of  Geophysical  Research. 


This  research  was  supported  by  the  National  Aeronautics  and  Spare 
Administration  through  grant  NGL-16-001-043  and  by  the  Office  of  Naval 
Research  through  grant  N00014-76-C-0016. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (Wf ,.n  Oaf*  Entered) 


REPORT  OOCUMENTATIOH  PAGE 


I.  REPORT  NUMBER 

U.  Of  1  IIW.I  84-9 


4.  TITLE  (and  Submit) 

I'ATII  I  NTF.C  RATED  CKOWT1I  OK  AURORAL  KILOMETR1C 
RADIATION 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


J  RECIPIENT'S  CATALOG  NUMBER 


S  TYPE  OF  REPORT  A  PERIOO  COVERED 

Progress  Februarv  1984 


6  PERFORMING  ORG.  REPORT  NUMBER 


7.  AuTHORfaJ 

N.  OM 1 1)1  and  I).  A.  Ol.’RNETT 


6.  CONTRACT  OR  GRANT  NUMBERfaJ 

N00014-76-C-0016 


».  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  '0.  PROGRAM  ELEMENT.  PROJECT,  TASK 

,  1,1  ■  a  AREA  *  WORK  UNIT  NUMBERS 

Depart  incut  ot  Physics  and  Astronomy 
University  of  Iowa 
Iowa  City,  I A  52242 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Electronics  Program  Office 
Office  of  Naval  Research 
Arlington,  VA  22217 


14  MONITORING  AGENCY  NAME  A  ADDRESS  (It  dltlerent  Item  Controlling  Oil  lea;  IS.  SECURITY  CLASS,  (el  thle  report) 


12.  REPORT  DATE 

February  1984 


IS.  NUMBER  OF  PAGES 

48 


UNCLASSIFIED 


ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  (ol  Ihl a  RaportJ 


A|)|>n>ved  lot  public  release;  distribution  is  unlimited. 


17  DISTRIBUTION  STATEMENT  (ol  lha  abalracl  entered  In  Block  30,  II  dlllermnt  (root  Report) 


l».  SUPPL  EMENTARY  NOTES 

Submitted  to  JOURNAL  OF  GEOPHYSICAL  RESEARCH 


19.  KEY  WORDS  (Continum  on  rovoroo  mid*  it  n*ca**arr  and  Identity  by  block  numb  or) 

Auroral  Kilometric  Radiation 
Ray  Tracing 


20.  ABSTRACT  (Contlnu*  on  rovono  of  do  it  n*c****ry  mnd  identify  by  btock  numbmr) 


(S<‘c  following  page) 


1  J  AN  7)  1473  COITION  OF  •  NOV  «>  It  OBSOLETE 
S/N  0102-LF-014-6601 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Afcon  Data 


ABSTRACT 


*4 1  : 

— ;•  Using  Poeverlein's  graphical  method,  three  dimensional  ray  path 
calculations  are  performed  to  evaluate  the  path-integrated  growth  of 
auroral  kilometric  radiation  (AKR).  The  ray  tracing  results  indicate 
that  waves  whose  initial  wave  vector  lie  in  the  local  meridian  plane 
continue  to  propagate  in  that  plane  and  that  among  these  waves ,  those 
with  frequencies  near  the  cutoff  frequency  (fg=o)  refract  substan¬ 
tially,  whereas  those  with  frequencies  well  above  the  cutoff  frequency 
suffer  little  refraction.  It  is  also  shown  that  waves  whose  initial 
wave  vector  lie  outside  of  the  local  meridian  plane  propagate  in  < he 
longitudinal  as  well  as  the  radial  and  the  latitudinal  directions.  The 
refraction  of  these  waves  is  also  highly  dependent  upon  the  wave  fre¬ 
quency,  i.e.,  waves  with  frequencies  near  fR=o  refract  substantially, 
whereas  waves  with  frequencies  much  above  ffc=o  undergo  little  refra< — 
tion.  In  order  to  test  the  electron  cyclotron  maser  mechanism  as  a 
method  for  generation  of  AKR,  a  typical  electron  distribution  function 
measured  in  the  auroral  zone  by  the  S3-3  satellite,  is  used  to  calcu¬ 
late  path-integrated  growths  of  representative  rays.  The  results  of 
this  study  indicate  that  electron  distribution  functions  like  those 
measured  by  the  S3-3  satellite  are  not  capable  of  amplifying  cosmic 


noise  background  to  the  observed  intensities  of  auroral  ki lorn-trie 


radiation,  and  that  much  steeper  slopes  { at  the  edges  of  the  loss 

i>V  i 
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cone  are  required.  The  presence  of  such  distribution  functions  in  the 
auroral  zone  is  plausible  If  one  assumes  that  backseat tered  electrons 
in  this  region  have  energies  less  than  a  few  hundred  eV. 


I.  INTRODUCTION 


Among  the  many  theories  proposed  for  the  generation  of  auroral 
kilometric  radiation  (AK.R) ,  the  one  suggested  by  Wu  and  bee  [1979]  has 
been  most  consistent  with  the  particle  and  wave  observations  in  the 
auroral  zone.  In  this  theory  mildly  relativistic  electrons  with  a  loss 
cone  distribution  amplify  right-hand  polarized  extraordinary  ( R— X ) 
waves  via  the  cyclotron  maser  mechanism.  Observations  of  auroral  kilo- 
metric  radiation  indicate  that  AKR  consists  of  both  right-  and  Icll- 
hand  polarized  radiation,  with  the  most  intense  and  dominant  component 
of  this  radiation  in  the  R-X  mode  [Gurnett  and  Green,  1978;  Kaiser  ct 
al.,  1978;  Benson  and  Calvert,  1979;  Shawhan  and  Gurnett,  1982). 
Measurements  of  the  electron  pitch  angle  distributions  in  the  auroral 
zone  by  the  S3-3  satellite  have  confirmed  the  existence  of  loss  cones 
in  the  upgoing  side  of  these  distributions  [Mizera  and  Fennell,  1977; 
Croley  et  al.,  1978]. 

In  the  past  few  years  much  work  has  been  done  to  further  develop 
the  work  done  by  Wu  and  Lee  [1979]  and  improve  our  understanding  of  Hie 
cyclotron  maser  mechanism  [Lee  and  Wu,  1980;  Omidl  and  Gurnett,  1987; 

Wu  et  al.,  1982;  Melrose  et  al.,  1982;  Dusenbery  and  Lyons,  1982; 

Hewitt  et  al.,  1982;  Hewitt  and  Melrose,  1983;  LeQueau  et  al . ,  1983; 
Omldi  et  at.,  1984] .  In  these  works  the  relativistic  cyclotron  reso¬ 
nance  condition  has  been  extensively  studied  and  its  properties  are  now 
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well  understood.  Furthermore,  in  order  to  improve  the  growth  rate 
calculations  done  by  Wu  and  Lee  [1979],  more  accurate  expressions  for 
the  growth  rate,  and  more  representative  distribution  functions  have 
been  used.  For  example,  Omidi  and  Gurnett  [1982]  and  Oraidi  et  al. 

[ 1984]  used  a  typical  electron  distribution  function  measured  by  the 
S3- 3  satellite  to  calculate  growth  rates  of  AKR.  Wu  et  al.  [1982]  and 
Melrose  et  al.  [1982]  have  also  computed  growth  rates  by  using  model 
electron  distribution  functions  based  on  measurements  of  the  S3-3  and 
ERSO  4  spacecrafts,  respectively.  The  results  of  these  growth  rate 
calculations  have  shown  that  waves  with  frequencies  near  the  electron 
gyrof requency  and  large  wave  normal  angles,  can  grow  to  large  ampli¬ 
tudes  by  extracting  the  free  energy  of  the  electrons  in  the  loss  cone. 

However,  as  was  pointed  out  by  Omidi  et  al.  [1984],  the  fact  that 
only  waves  with  large  wave  normal  angles  can  be  amplified  has  raised 
the  question  as  to  whether  refraction  of  waves  in  the  source  region  can 
hamper  their  growth.  In  order  to  answer  this  question,  it  is  necessary 
to  calculate  the  path  integrated  growth  of  these  waves  and  see  whether 
an  overall  amplification  is  obtained  or  not.  In  this  paper,  three 
dimensional  ray  path  calculations  of  the  auroral  kilonetric  radiation 
are  performed  by  using  Poeverlein's  graphical  method.  This  method 
employs  the  fact  that  as  rays  propagate  in  a  plain  stratified  medium, 
the  components  of  the  wave  vector  (k)  in  the  stratified  plane  remain 
constant.  The  results  of  this  ray  tracing  study  indicate  that  waves 
with  frequencies  (f)  near  the  local  R-X  cutoff  frequency  (fR»o)  suffer 


the  most  refraction  and  as  f  gets  larger  the  amount  of  refraction 


becomes  less. 

In  order  to  compute  the  path  Integrated  growth  of  the  waves,  it 
has  been  assumed  that  the  electron  pitch  angle  distribution  remains  the 
same  along  the  ray  path,  and  therefore  the  only  parameters  that  change 
are  the  electron  plasma  frequency  (fp),  the  electron  gyrof requency 
(fg),  the  group  velocity  (vg),  and  the  wave  normal  angle  (  j).  The 
expression  for  the  growth  rate  and  the  electron  distribution  function 
used  in  computing  the  growth  rates  are  the  same  as  those  given  in  Omidi 
et  ai.  |  198'* )  . 

Computations  of  the  path  integrated  growth  of  AKK  have  shown  Dial 

while  the  S3-3  electron  distribution  functions  are  unstable,  they  are 

not  capable  of  providing  sufficient  amplification  for  the  kilomotric 

radiation.  Using  two  modified  electron  distribution  functions  witli 

JF 

much  steeper  slopes  (-^j— )  at  the  edges  of  the  loss  cone,  it  is 
demonstrated  that  the  cyclotron  maser  mechanism  can  be  a  viable  genera¬ 
tion  mechanism  for  AKR  if  electron  distributions  with  larger  slopes 
DF 

(-r — )  exist  in  the  auroral  region.  Since  measurements  of  sharp  edges 
Jv  i 

in  the  distribution  functions  are  extremely  difficult,  one  cannot  rule 
out  the  possibility  of  their  existence  based  on  the  present  spacecraft 
data.  Some  elementary  theoretical  considerations  seem  to  indicate  that 
electron  distribution  functions  with  much  steeper  slopes  at  the  edge  of 
the  loss  cone  may  exist  if  one  assumes  that  all  the  backseat tered  elec¬ 
trons  have  energies  less  than  a  few  hundred  eVs. 
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II.  RAY  TRACING 

In  this  section  Poeverlein's  graphic  method  of  ray  tracing  is 
briefly  discussed,  and  the  results  of  the  three-dimensional  ray 
tracings  for  the  auroral  kiloraetric  radiation  are  presented. 

In  Poeverlein's  graphic  method  of  ray  tracing  it  is  assumed  that 

waves  propagate  in  a  plane  stratified  medium.  This  method  utilizes  the 

fact  that  although  the  wave  vector  k  changes  along  the  ray  path,  the 
► 

components  of  k  in  the  stratified  plane  remain  constant. 

For  the  model  used  in  this  paper,  it  is  assumed  that  the  stratified 
plane  (x,  y)  is  perpendicular  to  the  gradient  of  the  magnetic  field 
(7B)  as  shown  in  Figure  la.  This  model  is  justified  by  noting  thal 
since  in  the  auroral  region  the  fp/fg  ratio  is  usually  very  small 
(fp/fg  <  .1),  [Benson  and  Calvert,  1979],  small  variations  of  the 
electron  density  can  he  ignored.  The  refraction  of  the  waves  is  then 
mainly  determined  by  the  geomagnetic  field,  which  is  assumed  to  be 
dipolar.  In  Figure  la,  (i,  j,  k)  forms  a  Cartesian  coordinate  system 
in  which  k  is  along  the  magnetic  pole  and  (i,  j)  lie  in  the  magnetic 
equator.  The  origin  of  the  (x ,  y,  z)  system  is  the  launching  point  for 
all  the  rays  and  in  this  study  it  is  located  at  a  radial  distance  of 
R  =  13,940  km,  a  magnetic  latitude  of  A  =  70°  and  a  magnetic  longitude 
of  1,  =  0°.  Note  that  the  (y,  z)  plane  lies  on  the  (k,  f )  plane  which 
del  lues  a  magnetic  meridian. 
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Since  the  refraction  of  waves  is  assumed  to  be  controlled  by  the 

Karlh's  dipole  magnetic  field,  the  model  for  the  electron  density  used 

In  the  present  study  is  a  simple  one.  The  fp/fg  ratio  is  kept  at  0.05 

throughout  the  region  where  the  ray  tracing  is  performed.  This  model 

Ik  consistent  with  the  fact  that  gradients  In  the  electron  densltv 

should  have  a  negligible  effect  on  the  ray  path  when  fp/fg  ratio  is 

small.  It  should,  however,  be  mentioned  that  this  simple  model  ignores 

the  density  enhancements  within  the  auroral  plasma  cavity  reported  by 

Benson  and  Calvert  [1979].  These  density  enhancements  could  result  in 

partial  reflection  of  the  waves,  as  suggested  by  Calvert  [1982]  in  his 

laser  feedback  model  of  AKR.  We  will  come  back  to  this  point  later. 

As  mentioned  earlier,  the  origin  of  the  (x,  v,  z)  system  is  the 

launching  point  for  all  rays.  Once  the  wave  frequency  and  the 
► 

direction  of  k  are  chosen,  the  components  of  the  Index  of  refraction 

y  l'  k.  '*  ** 

vector,  N  =  —  along  the  x  axis  (Nxo)  and  the  y  axis  (Nyo)  are 
determined  and  remain  unchanged  along  the  ray  path.  As  the  ray  is 
stepped  forward  along  the  group  velocity  (which  is  perpendicular  to 
the  index  of  refraction  surface)  a  new  value  for  the  index  of  refrac¬ 
tion  is  computed.  By  constructing  a  new  coordinate  system  (x',  y',  z') 
parallel  to  (x,  y,  z),  with  its  origin  located  at  the  new  position  of 
the  ray,  new  values  of  the  wave  normal  angle  ii  and  the  azimuthal  angle 
r  (see  Figure  lb)  can  be  determined  such  that  the  components  of  N  in 
the  x'  and  y'  direction  are  Nxo  and  Ny0,  respectively.  Note  that 

♦  “  A 

although  at  the  Launching  point,  B  is  in  the  (y,  z)  plane,  at 

> 

subsequent,  points  along  the  ray  path  the  local  B  will  not  lie  in  the 
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(y',  z')  plane  unless  initially  =  0,  in  which  case  the  ray  propagates 
in  the  meridian  plane. 

Before  illustrating  specific  rays  that  were  computed  in  this 
study,  it  is  helpful  to  review  the  general  results  which  are  summarized 
in  Figure  2.  This  figure  shows  the  ray  paths  for  four  waves  with 
different  frequencies  and  initial  k  direction,  all  starting  from  the 
same  point.  Rays  (a)  and  (b)  lie  in  the  front  meridian  plane  (i.e., 
initial  =  0  and  k  has  no  longitudinal  component).  Ray  (a)  lias  a  fre¬ 
quency  near  the  cutoff  frequency  and  is  refracted  considerably,  whereas 
ray  (b)  which  lias  a  frequency  much  above  the  cutoff  frequency  sutlers 
much  less  refraction.  This  dependence  is  easily  explained  by  noting, 
that  initially,  the  index  of  refraction  N  of  ray  (a)  is  very  small  -(due 
to  the  fact  that  f  is  near  f^-g)  anc*  therefore  the  components  of  N  in 
the  stratified  plane  are  very  small.  As  ray  (a)  propagates,  N  gets 
larger  and  in  order  for  the  components  of  N  in  the  stratified  plane  to 

remain  constant,  the  wave  normal  angle  -3  has  to  decrease  (note  that 
■+1 

B  is  nearly  perpendicular  to  the  stratified  plane,  i.e.,  ^  in  Figure  lb 
is  small).  Ray  (b)  on  the  other  hand  has  N  -  1  and  its  index  of 
refraction  does  not  change  appreciably  as  it  propagates.  Therefore,  it 
will  not  suffer  as  much  refraction.  Rays  (c)  and  (d)  whose  Initial  k 
do  not  lie  in  the  front  meridian  plane  (initially  =  90°)  propagate  In 
three  dimensions  from  the  front  towards  the  back  meridian  plane.  Again 
ray  (c)  with  a  frequency  near  the  R-X  cutoff  frequency  refracts  sub¬ 
stantially,  whereas  ray  (d)  which  has  a  frequency  well  above  the  cutoff 
frequency  travels  longitudinally  without  too  much  refraction.  As  can 
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0c  soon  in  Figure  2,  rays  (a)  and  (b)  which  Initially  have 
propagate  only  in  the  radial  and  latitudinal  directions,  whereas  rnvs 
(c)  and  (d)  propagate  mostly  in  the  longitudinal  and  radial  directions. 
Thus,  as  the  initial  value  of  ,  increases,  the  rays  travel  more  in  the 
longitudinal  direction  and  less  in  the  latitudinal  direction. 

Having  described  the  general  behavior  of  the  R-X  node  waves  in  the 
AKR  source  region,  we  now  focus  on  six  specific  waves  whose  frequency 
and  initial  k  direction  are  shown  in  Table  1.  The  ray  paths  of  these 
waves  are  representative  of  the  different  types  of  ray  paths  that  are 
observed.  In  Figure  3  the  variations  of  the  wave  normal  angle  u  as  a 
function  of  path  length  are  shown.  The  total  path  length  for  each  ray 
shown  in  this  figure  is  100  km.  From  Figure  3  one  can  draw  the  follow¬ 
ing  conclusions.  Clearly  rays  #1  and  #2  whose  frequencies  are  near  the 
cutoff  frequency  fg-Q ,  suffer  substantial  refraction,  whereas  rays  'f3 
through  #6  whose  frequencies  are  well  above  fg=o>  do  not:  suffer  as  much 
refraction.  It  Is  also  evident  that  much  of  the  refraction  occurs  in 
the  first  twenty  kilometers  of  the  ray  path.  As  will  be  seen  later, 
this  effect  plays  a  crucial  role  in  preventing  waves  from  getting 
amplified.  Another  point  to  notice  in  Figure  3  is  that  two  rays  with 
the  same  frequency  but  different  initial  azimuthal  directions  <  are 
refracted  differently,  with  the  wave  propagating  in  the  meridian  plane 
suffering  more  refraction.  The  reason  is  that  for  a  ray  travelling  in 
the  meridian  plane  the  magnitude  of  B  changes  more  rapidly  than  for  a 
ray  propagating  in  the  longitudinal  direction.  However,  this  differ¬ 
ence  in  refraction  is  not  substantial,  and  as  will  be  shown  later  it 


does  not  lead  to  markedly  different  path  integrated  growths. 


In  Figure  4  the  variations  of  the  magnetic  longitude  and  latitude 
of  each  of  the  six  rays  for  a  path  length  of  400  km  are  shown.  The 
starting  point  for  all  the  rays  is  at  a  latitude  of  70°  and  a  long.it  ud 
of  0°.  The  fact  that  rays  // 2  and  #4  are  propagating  in  the  meridian 
plane  is  obvious.  Note  that  since  ray  f/2  undergoes  more  refraction 
than  ray  #4,  it  does  not  extend  to  as  low  a  latitude  as  does  ray  #4. 

In  other  words,  ray  #2  refracts  upward  much  more  than  ray  #4.  It  is 
evident  in  Figure  4,  that  rays  #5,  #3,  and  #1  with  their  initial  ;  = 
90°,  travel  mostly  in  the  longitudinal  direction  with  ray  # 1  propagat¬ 
ing  a  bit  more  in  the  latitudinal  direction  due  to  more  refraction. 
Finally,  ray  #6  whose  initial  azimuthal  angle  is  =  20°,  travels  both 
in  the  longitudinal  as  well  as  the  latitudinal  directions. 

Having  seen  the  general  behavior  of  the  waves  in  the  auroral 
region,  we  now  proceed  to  evaluate  the  growth  of  these  waves  by  the 
cyclotron  maser  mechanism. 


III.  PATH  INTEGRATED  GROWTHS  USING  MEASURED  DISTRIBUTION  FUNCTIONS 


As  was  mentioned  earlier,  although  the  electron  distribution  func¬ 
tions  measured  in  the  auroral  region  are  unstable  and  have  large  growth 
rates  in  certain  frequency  ranges,  a  large  growth  rate  by  itself  is  not. 
sufficient  to  insure  that  the  waves  will  grow  to  sufficiently  large 
amplitudes  along  the  ray  path.  To  demonstrate  this  point,  the  growth 
and  damping  rates  of  AK.R  as  a  function  of  frequency  for  four  different 
wave  normal  angles  are  shown  in  Figure  b.  These  growth  rates  were 
computed  by  using,  the  eLectron  distribution  function  measured  by  the 
S3-3  satellite  (see  Figure  6)  and  the  method  of  calculation  given  in 
Otaidi  et  al.  [1983],  From  Figure  5  it  is  clear  that  substantial  growth 
rates  can  be  obtained  only  when  b  ~  80°.  When  0  <  70°  all  waves  are 

damped.  The  reason  is  that  as  G  gets  smaller  most  of  the  resonant 

DF 

electrons  lie  outside  the  loss  cone  region  where  (— — )  is  negative 

1 

(see  Omidi  and  Gurnett  [1982]).  Thus,  the  only  way  that  waves  can  grow 
to  large  amplitudes  is  if  the  ray  paths  are  such  that  )  remains  close 
to  80°  for  a  considerable  length  of  time.  To  see  if  such  waves  indeed 
occur,  we  have  calculated  growth  factor  (G)  defined  as: 


G  exp [g] 


(1) 


»  - 
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whore  g  -  j 
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,L 


v  (■■:') 
g 


d  C ' 


Is  the  exponential  gain  factor.  (  1 '  )  and  Vg(..')  are  the  imaginary 
part  of  the  frequency  and  the  group  velocity  as  a  function  of  distance, 
respectively,  and  L  is  the  path  length.  The  growth  rates  jj(i')  are 
computed  in  the  same  manner  as  in  Omidi  et  al.  [1984],  The  group 
velocity  v^  is  given  by 


v 
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where  (~^)  and  (“)  are  found  by  using  the  dispersion  relation  obtained 
in  the  cold  plasma  approximation.  The  value  of  Vg  is  zero  when 
f  =  fg=o  and  rapidly  approaches  the  speed  of  light  as  If  :  (f  -  fR_o) 
increases  (see  Dusenbery  and  Lyons  [1982)). 

Intensity  measurements  of  AKR  indicate  that  the  growtli  of  the 
cosmic  noise  background  electric  field  by  a  factor  of  G  =  e  1  lj  Is 
sufficient  to  account  for  the  observed  intensities  of  the  kilonetric 
radiation,  i.e.,  g  ~  10  gives  the  required  amplification.  In  Figure  7, 
a  plot  g  vs.  L  for  the  six  different  rays  discussed  in  the  last  section 
are  shown.  In  order  to  calculate  g,  It  has  been  assumed  that  electron 
distribution  functions  similar  to  that  shown  In  Figure  6,  exist  along 
the  ray  path.  Figure  7  clearly  illustrates  that  rays  #1  and  #2  with 
frequencies  near  the  cutoff,  suffer  heavy  damping  due  to  the  rapid 
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refraction  that  they  undergo.  In  other  words,  although  these  rivs 
start  with  0  sufficiently  near  90°  for  large  growth  rates  they  i annul 
maintain  these  wave  normal  angles  long  enough  to  be  amplified.  Also 
note  that  although  ray  111  propagates  in  the  longitudinal  direction,  its 
damping  is  similar  to  that  of  ray  #2  which  propagates  in  the  meridian 
plane . 

As  can  be  seen  in  Figure  7,  rays  #3  and  *14  which  suffer  much  less 
refraction  are  InitialLy  amplified  but  are  later  absorbed.  Tlu-  Maximum 
g  for  these  rays  is  less  than  0.5.  Therefore,  the  initial  growths  of 
these  rays  are  not  anywhere  close  to  the  growth  required  to  account  I  or¬ 
tho  intensities  of  AKR.'  Similarly,  rays  lf 5  and  #6  are  Initially  ampli¬ 
fied  and  suffer  subsequent  damping.  These  rays,  however,  are  amplified 
less  and  propagate  a  larger  distance  before  getting  damped.  A  compar¬ 
ison  of  g  between  rays  It 3  and  It 4  or  It 5  and  itb  demonstrates  that  rays 
propagating  in  the  longitudinal  direction  travel  a  larger  distance 
before  getting  absorbed.  This  occurs  because  for  a  given  frequency 
waves  propagating  in  the  longitudinal  direction  suffer  less  refraction. 
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IV.  PATH  INTEGRATED  GROWTHS  WITH  ARTIFICIALLY 
STEEPENED  VELOCITY  GRADIENTS 

Judging  from  the  results  shown  in  Figure  7,  it  seems  certain  that 
with  electron  distribution  functions  like  those  measured  by  the  SI-3 
satellite,  the  cyclotron  maser  mechanism  cannot  account  for  the 
generation  of  AKR.  Thus,  if  the  cyclotron  maser  Instability  is  to 
remain  a  viable  mechanism  for  the  generation  of  AKR  it  is  necessary  for 
the  electron  distribution  in  the  loss  cone  to  have  steeper  slopes  (—■•) 
than  those  measured  by  S3-3  satellite.  Since  wave  particle  inter¬ 
actions  tend  to  fill  in  the  loss  cone  via  pitch  angle  scattering,  thus 
reducing  the  velocity  space  gradients,  it  is  possible  that  electron 
distribution  functions  with  steeper  slopes  do  in  fact  exist.  Further¬ 
more,  because  of  the  limitations  imposed  by  the  angular  (fields  of  view 
are  7°  *  10°  FWHfl)  and  temporal  (20  sec.  for  a  complete  angular  distri¬ 
bution)  resolution  of  the  S3-3  particle  instrumentation,  it  is  possible 
that  electron  distribution  functions  with  velocity  gradients  much 
steeper  than  those  shown  in  Figure  6  exist,  but  are  not  detected.  In 
this  section  we  estimate  the  minimum  velocity  space  gradients  (n  the 
loss  cone  that  would  be  required  to  amplify  the  cosmic  noise  background 
to  the  observed  intensities  of  AKR. 
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To  construct  electron  distribution  functions  with  steeper  slopes 
(-7^--)  in  the  loss  cone  than  those  measured  bv  the  S3-3  satellite,  we 
used  the  following  approach.  Generally  speaking,  the  electrons  In  the 
loss  cone  can  he  broken  into  three  separate  populations,  one  population 
consisting  of  the  electrons  that  were  pitch  angle  scattered  into  the 
loss  cone  via  the  wave  particle  interact  tons ;  the  second  population 
bei ng  composed  of  the  backseat tered  electrons;  and  finally,  the  third 
population  consisting  of  the  electrons  with  ionospheric  origin.  The 
electron  distribution  function  shown  in  Figure  6  was  measured  in  the 
auroral  region  and  at  the  radial  distance  of  about  2.2  Rp  where  an 
electrostatic  potential  drop  of  about  1  kV  with  respect  to  the  top  of 
the  Ionosphere  Is  believed  to  exist  [Croley  et  al .  ,  1978).  Since  most 
ot  the  Ionospheric  electrons  have  energies  less  than  1  keV,  then  one 
can  assume  that  these  electrons  will  not  be  able  to  reach  the  S 1—  3 
satellite.  As  for  the  backscattered  electrons,  a  detailed  study  of  the 
collisions  of  the  downgoing  electrons  with  the  particles  In  the  iono¬ 
sphere  is  required  to  determine  their  energy  spectrum.  Since  such  a 
study  is  beyond  the  scope  of  this  paper,  we  simply  assume  that  the 
backscattered  electrons  at  the  height  of  the  S3-3  satellite  have  ener¬ 
gies  of  a  few  hundred  eVs  or  less  so  that  they  do  not  resonate  with 

waves.  Based  on  this  assumption,  all  the  high  energy  electrons  in  the 

loss  cone  are  due  to  pitch  angle  scattering.  Therefore,  the  electron 
dlst rlbut Ion  function  prior  to  wave  particle  interactions  must  have  had 
an  empty  loss  cone.  However,  even  if  tne  loss  cone  is  empty,  a  limita¬ 
tion  still  exists  on  the  steepness  of  the  loss  cone  because  the 


absorption  in  the  atmosphere  is  not  an  abrupt  boundary.  This  limita¬ 


tion  may  be  analyzed  in  the  following  way.  At  high  altitudes  the  elec¬ 
trons  are  collisionless,  however,  as  they  approach  the  ionosphere  they 
reach  a  point  on  the  field  line  (Mj)  where  collisions  with  particles  in 
the  atmosphere  start  to  become  important.  This  point  can  be  taken  to 
be  the  mirroring  point  for  the  electrons  and  defines  the  loss  cone 
angle  aj .  It  is  obvious,  however,  that  not  all  the  electrons  passing 

this  point  wilL  immediately  suffer  a  collision,  and  therefore  some  of 

these  electrons  will  mirror  back  before  suffering  a  collision.  Wo  now 
define  a  second  point  on  the  same  field  line  and  below  M|  such  that 
all  the  electrons  passing  Mj  wilt  either  he  mirrored  hack  or  scattered 
before  reaching  1^,  i.e.,  the  distance  between  Mj  and  M2  is  on  the 
order  of  an  electron  mean  free  path.  We  denote  this  distance  by  AM. 
Point  M2  defines  a  loss  cone  angle  02  below  which  essentially  no  par¬ 
ticles  exist  in  the  loss  cone.  Since  some  electrons  passing  the  point 
Mj  are  scattered  while  others  are  i rrored  before  reaching  M2,  one  can 
see  that  the  electron  distribution  in  the  loss  cone  drops  abruptly  as 
one  goes  from  aj  to  u2*  This  scattering  process  can  be  used  to  put  a 
limitation  on  the  steepness  of  the  slopes  at  the  edge  of  the  loss  cone. 
Clearly  as  Mj  and  M2  approach  each  other  then  )  becomes  Infinite 

and  as  the  distance  between  Mj  and  M2  increases  the  gradients  become 

less  steep. 

To  properly  construct  a  model  electron  distribution  function,  it 


Is  necessary  to  consider  the  effect  of  a  parallel  electric  field. 
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Persson  [1966]  user!  the  conservation  of  energy  and  the  first 
odi.ih.il  i  c  Invariant  to  show  that  in  tin1  presence  ot  parallel  electric 
fields,  the  loss  cone  boundary  is  defined  by  a  hyperbola  (see  Figure  6) 
g i veu  by 


H 

'¥ 


V)  =  0 


(3) 


where  v h  and  vt  are  the  electron  velocity  components  parallel  and  per¬ 
pendicular  to  the  magnetic  field,  respectively;  B;  and  B  are  the  mag¬ 
netic  field  at  mirror  point  and  the  observation  point;  V;.  and  V  are  the 
electric  potential  at  the  mirror  and  the  observation  points,  respec¬ 
tively;  in  is  the  electron  mass,  and  e  is  its  charge.  Since  B  is 
measured  by  the  S3-3  satellite,  one  only  needs  to  know  B  ,•  and  .*  V  (  v 

-V')  to  do l ermine  the  loss  cone  boundary  given  by  Equation  1.  As  was 
mentioned  earlier,  typically  AV  =  IkV.  Location  of  Mj  on  the  other 
hand  is  not  exactly  known  and  it  can  vary  between  I  no  to  300  km  above 
the  Earth's  surface.  Note,  however,  that  since  By  does  not  change 
appreciably  within  this  height  interval,  the  loss  cone  boundary  corre¬ 
sponding  to  Mj  can  be  adequately  determined  by  assuming  ?lj  to  be  at  an 
altitude  of  300  km.  This  boundary  is  shown  in  Figure  6. 

In  order  to  construct  an  electron  distribution  function  F( v ,  a) 
(where  v  is  electron  velocity)  with  steeper  slopes  in  the  loss  cone 
than  those  In  the  S3-3  distribution  function,  point  M2  is  chosen  and 
the  loss  cone  boundary  corresponding  to  this  point  is  determined.  It 
Is  then  assumed  that  F(v,u)  is  zero  on  and  below  the  loss  cone  boundary 


of  M2  and  that  it  rises  linearly  with  v  and  a  to  the  values  of  the  Si—  i 
distribution  on  the  loss  cone  boundary  of  Mj .  F(v,u)  outside  the  loss 
cone  is  the  same  as  the  S3-3  distribution  function.  The  assumption 
that  F(v,u)  rises  linearly  with  v  and  a  between  the  two  loss  cone 
boundaries  can  be  justified  by  noting  that  these  two  boundaries  are 
extremely  close  to  each  other  and,  thus,  in  an  expansion  of  K(v,0 
between  these  two  boundaries  the  second  and  higher  order  terms  can  be 
ignored. 

In  order  to  find  the  minimum  slopes  in  the  loss  cone  that  are 
required  to  result  in  a  gain  factor  of  g  ~  10,  AM  has  been  varied 
between  3  and  100  km.  An  estimate  of  a  reasonable  value  of  AM  wi 1 1  be 
given  later.  In  Figure  8  the  growth  rates  resulting  from  a  distribu¬ 
tion  function  with  AM  =  3  km  are  shown  (run  I* 2).  Thu  difference 
between  these  growth  rates  and  those  shown  in  Figure  S  is  striking.  As 
can  be  seen,  large  growtli  rates  are  obtained  even  when  u  <  70°.  The 

reason  is  that,  although  for  J  <  70°  most  of  the  resonating  electrons 

c)F 

are  outside  the  loss  cone  region  where  (-j^— )  is  negative,  the  large  and 

IF  1 

positive  (~~)  at  the  edge  of  the  loss  cone  can  easily  overcome  these 
negative  contributions,  and  result  in  an  overall  growtli  for  the  waves. 
Note  that  u>^  =  0  when  the  resonance  ellipse  falls  below  the  loss  cone 
boundary  corresponding  to  point  M2,  because  there  are  no  electrons 
below  this  boundary.  Note  also  that  the  two  loss  cone  boundaries 
corresponding  to  Mj  and  M2  are  so  close  that  they  could  not  he  resolved 
in  Figure  b,  and  that  at  a  given  velocity  the  difference  In  the  loss 
cone  angles  A  a  =  uj  -  02  is  on  the  order  of  10-^  degrees. 


Tin-  pat  1)  i  nt  egrat  ed  growths  ot  rays  /M  t  ii roujOi  #h  for  t  ho  case  of 
AM  =  3  km  (run  It2)  are  shown  in  Figure  H.  It  ran  ho  sivn  that  rays  ?  1 
and  l>2  whoso  group  velocity  Is  much  loss  than  the  speed  of  light  are 
greatly  amplified  within  a  very  short  distance  (L  -  A  km).  On  the 
other  hand,  rays  #3  through  It b  with  v^,  closer  to  the  speed  of  light 
take  a  longer  distance  for  a  similar  amplification.  At  any  rate,  it  is 
obvious  that  all  six  rays  are  sufficiently  amplified  (g  10)  within 
2  -  20  km,  and  thus  the  electron  distribution  function  modeled  by 
taking  .AM  =  3  km  is  easiLy  capable  of  generating  AKK. 

Similar  calculations  with  Mn  at  different  positions  have  shown 
that  the  distribution  funet Ion  with  the  minimum  required  slopes  In  t hi* 
loss  cone  Is  obtained  when  ‘11  =  /()  km  (run  /*3)  which  corresponds  to 
A»  •  0.S  degrees.  The  .  cowth  rates  for  this  rase  are  shown  In  Figures 
10  and  II  for  >  varying  from  80°  down  to  33''.  Note  that  In  this  case 
both  growth  and  damping  can  occvir,  depending  upon  the  wave  frequency. 

In  the  paper  by  Omidi  and  Gurnet t  |19K2]  a  qualitative  picture  was 
given  as  to  why  two  regions  of  growth  and  one  region  of  damping  are 
obtained  (see  Figure  3a),  when  the  S3-3  distribution  function  is  used 
to  calculate  .  It  was  show  i  that  for  wave  frequencies  within  the 
first  or  second  growth  regions  most  or  all  of  the  resonance  ellipses 
fall  inside  of  the  loss  cone,  whereas  at  frequencies  within  the  damping 
region,  most  of  the  resonance  ellipses  fall  outside  of  the  loss  cone, 
and  hence  waves  at  these  frequencies  are  damped.  However,  in  the  case 
ot  the  growth  rates  shown  in  Figures  10  and  II,  a  similar  argument  does 


not  applv  because  of  the  absence  of  any  electrons  below  the  loss  cone 
boundary  of  Mi ,  While  in  this  case  a  qualitative  picture  for  the 
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behavior  of  the  growth  rates  does  not  exist,  close  examination  of  the 
growth  rate  calculations  indicate  that  the  presence  of  trapped  elec¬ 
trons  in  the  distribution  function  plays  an  important  role  in  the  over 
all  growtii  rates  of  the  waves.  In  other  words,  while  (-—— )  is  large 
and  positive  at  the  edge  of  the  loss  cone,  it  is  not  large  enough  to 

always  overcome  the  negative  (— — )  outside  the  loss  cone,  and  thus 

3V1 

whether  the  resonance  ellipse  for  a  particular  frequency  falls  in  the 

)F 

regions  of  positive  or  negative  (— — )  in  the  trapped  part  of  the 
distribution  function,  can  affect  the  outcome  of  the  growth  rale 
calculations  for  that  frequency. 

The  path  integrated  growths  of  rays  #1  through  #6  correspond¬ 
ing  to  run  It3  are  shown  in  Figure  12.  In  this  case  rays  It  1  and  0 2  can 
be  sufficiently  amplified  (g  ~  10)  within  70  ktn.  but  on  the  other  hand 
rays  It3  and  #4  are  initially  amplified  and  suffer  subsequent  damping. 
The  reason  for  this  difference  can  be  seen  by  noting  that  for  large  Y 
fg/f  (Y  >  .99)  growth  can  only  be  obtained  when  0  -  80°  or  ‘j  <  80°. 
Since  rays  It  1  and  It 2  suffer  extensive  refraction  their  wave  normal 
angles  soon  become  less  than  60°  and  can  therefore  be  amplified.  The 
fact  that  these  rays  have  group  velocities  much  less  than  r  also  helps 
to  reduce  the  path  lengths  required  to  obtain  sufficient  ampl i f icar i on 
Rays  It 3  and  #4  on  the  other  hand  do  not  suffer  much  refraction  and 
therefore  their  wave  normal  angle  will  not  go  below  60°.  These  rays 
are  therefore  Initially  amplified  (when  L;  -  80°)  but  as  they  propagate 
they  begin  to  damp.  Rays  // 5  and  It6  suffer  even  less  refraction  than 
rays  It 3  and  #4.  Therefore,  It  takes  a  longer  distance  before  they 
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begin  to  damp.  It  is  clear  from  Figures  10  and  11  that  rays  with 
Y  <  .99  and  initial  0o  ~  75°  (not  traced  in  this  study)  may  also  be 
amplified.  However,  since  the  group  velocities  of  these  waves  would  he 
close  to  the  speed  of  light,  the  path  lengths  required  for  sufficient 
amplification  would  he  on  the  order  of  1000  km  which  is  much  too  long. 

As  was  mentioned  earlier,  the  smallest  velocity  gradients  in  the 
loss  cone  which  allow  for  stiff icient  amplification  of  AKR  correspond  to 

0.1  degrees.  To  see  if  such  a  .',.i  is  reasonable,  the  mean  free 
paths  of  the  electrons  ( ?,)  as  a  function  of  height  and  electron  energy 
have  been  computed.  These  calculations  were  performed  by  using 
ionization  cross  sections  of  O2 ,  ^2,  and  0  given  by  Equation  5  in  Banks 
et  al.  [1974]  at  electron  energies  of  0.5  and  5  keV.  The  atmospheric 
neutral  densities  were  obtained  from  the  COSPAR  International  Reference 
Atmosphere,  1972.  In  Figure  13,  two  plots  of  vs.  altitude  are  shown 
for  an  electron  energy  of  5  keV.  The  mean  free  paths  of  0.5  keV 
electrons  are  smaller  by  about  a  factor  of  5.  In  panel  (a)  of  Figure 
13,  the  "maximal"  model  of  the  atmosphere  was  used  to  determine  the 
neutral  densities  as  a  function  of  height,  whereas  in  panel  (b)  the 
"minimal"  model  was  utilized. 

Using  the  electron  mean-free  paths  shown  in  panel  (a),  we  have 
calculated  the  electron  reflection  coefficient  r  defined  as: 
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where  nr  is  the  density  of  the  reflected  electrons  at  a  given  energy 
and  pitch  angle  and  is  given  by: 


n 
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,  s  ds' 

V*p^2  i 


(■>) 


In  Equation  5,  nQ  is  the  electron  density  at  a  given  energy  and  pitch 
angle  at  the  height  where  the  electron  mean  free  path  is  essentially 
infinite  (in  this  study  altitude  of  700  km).  Also,  s  corresponds  to 
the  distance  travelled  by  the  electrons  along  their  helical  orbit 
before  they  are  mirrored  back.  The  factor  of  two  in  the  exponential 
accounts  for  the  scattering  of  the  electrons  after  they  are  mirrored 
back.  Note  that  because  the  electron  orbit  is  a  function  of  both 
energy  and  pitch  angle  and  also  is  a  fund  ion  of  energy,  r  is  a 
function  of  both  the  electron  energy  and  pitch  angle. 

In  Figure  14,  two  plots  of  r  vs.  pitch  angle  are  shown.  In  panel 
(a)  of  this  figure,  the  reflection  coefficients  of  0.5  keV  electrons 
are  shown,  whereas  in  panel  (b)  those  of  5  keV  electrons  are 
illustrated.  As  can  be  seen  in  Figure  14,  the  reflection  coefficient  r 
drops  rapidly  from  one  to  zero  within  a  couple  of  degrees  in  both 
panels  (a)  and  (b).  However,  the  pitch  angle  at  which  r  begins  to  drop 
is  a  function  of  electron  energy.  This  effect  is  due  to  the  parallel 
electric  field,  and  is  consistent  with  the  fact  that  in  the  presence  <>( 
a  parallel  electric  field,  the  loss  cone  boundary  Is  a  hyperbola.  At 
any  rale,  it  Is  clear  that,  using  the  electron  mean  free  paths  shown  In 
panel  (a)  of  Figure  13,  An  is  about  two  degrees.  However,  decreasing 
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tlx'  mean  f rt-o  paths  by  a  facto r  of  three  or  four  would  allow  r  to  drop 
from  ().*)  to  0.1  within  half  degrees  or  less.  Since  in  this  study  we 
have  only  considered  ionization  cross  sections,  and  have  ignored  other 
types  of  collisions  like  elastic  and  Coulomb  collisions,  it  is  possible 
that  electron  mean  free  paths  are  indeed  smaller  than  those  shown  in 
Figure  13.  In  that  case,  in  ~  0.3  degrees  could  be  a  reasonable  value. 


25 


V.  DISCUSSION  AND  CONCLUSION 

The  aim  of  this  study  has  been  to  see  if  a  typical  electron  dis¬ 
tribution  function  measured  in  the  auroral  region  by  the  S3-3  can 
sufficiently  amplify  R-X  mode  waves.  It  has  been  demonstrated  that 
the  electron  distribution  function  measured  by  the  S3-3  satellite 
(shown  in  Figure  6)  cannot  result  in  sufficient  amplification  of  waves, 
and  that  much  steeper  slopes  at  the  edges  of  the  loss  cone  are 
required.  By  artificially  steepening  the  velocity  space  gradients  at 
the  edge  of  the  loss  cone,  we  have  shown  that  the  angular  width  of  the 
edge  of  the  loss  cone  must  be  au  ~  0.5  degrees  in  order  to  have  ade¬ 
quate  growth.  Assuming  that  hackscattered  electrons  and  ionospheric 
electrons  are  negligible,  computations  of  the  absorption  of  electrons 
by  collisions  with  the  neutral  atmosphere  shows  that  the  edge  of  the 
loss  cone  can  have  a  width  as  small  as  Au  -  2  degrees.  However,  inclu¬ 
sion  of  the  elastic  and  Coulomb  collisions  in  the  electron  mean  free 
path  (~)  calculations  and  also  the  uncertainty  in  the  atmospheric 
neutral  densities  could  decrease  £  by  a  factor  of  3  or  more  which  would 
then  reduce  Au  to  about  0.5  degrees.  Thus,  loss  cone  edges  with  widths 
as  small  as  Au  ~  0.5  degrees  are  in  principle  possible.  Future 
observations  and  more  detailed  theoretical  study  of  the  backseat t eri ng 
of  downgoing  electrons  with  the  Ionospheric,  particles  are  needed  to 


establish  whether  such  steep  velocity  gradients  actually  exist.  At 
present,  no  measurements  are  available  with  sufficiently  good  angular 


reso lut  ion  to  resolve  such  steep  gradients.  Measurements  of  this  type 
are  needed  and  should  he  made  on  future  auroral  missions. 

finally,  a  few  comments  on  the  laser  feedback  model  proposed  for 
i  he  gene rat  ion  of  AKK  by  Calvert  (1982)  are  in  order.  In  this  model 
density  Irregularities  are  employed  to  partially  reflect  the  waves  and 
thereby  allow  an  oscillator  type  mechanism  to  operate.  It  has  been 
pointed  out  by  Calvert  ( 1982]  that  in  this  model,  a  total  wave  growth 
of  about  40  db  per  loop  is  required.  It  is  clear  from  Figure  7  that 
again  such  gains  cannot  be  achieved  by  the  electron  distribution  func¬ 
tion  shown  in  Figure  6.  In  addition,  since  a  wave  propagating  along  a 
loop  will  at  some  portions  of  the  loop  (downgoing  side)  resonate  with 
the  downgoing  electrons  and  most  probably  be  damped  (see  Omidi  et  al., 
|1484]),  the  amplification  along  the  upgoing  portion  of  the  loop  will 
have  to  compensate  for  this  damping.  This  again  demonstrates  the  need 
for  steeper  slopes  In  the  loss  cone,  if  the  cyclotron  maser  mechanism 
is  to  remain  a  viable  one. 
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Figure  1 


Figure  2 


Figure  3 


FIGURE  CAPTIONS 

In  la,  the  coordinate  systems  (i,  j,  k)  with  k  along 
the  magnetic  pole  and  (x,  y,  z)  whose  origin  is  tin* 
launching  point  for  all  rays  are  shown.  Note  that 
(x,  y)  plane  forms  the  stratified  plane.  In  Ih  the 
index  of  refraction  surface  and  the  wave  normal  angle 
j  and  the  azimuthal  angle  are  shown.  (x‘,  y',  z') 
is  parallel  to  (x,  y,  z). 

Typical  ray  paths  for  four  waves  are  shown.  Kays  (a) 
and  (b)  propagate  in  the  front  meridian  plane, 
whereas  rays  (c)  and  (d)  travel  in  the  longitudinal 
as  well  as  radial  and  latitudinal  directions.  Rays 
(a)  and  (c)  have  frequencies  near  the  cutoff  whereas 
rays  (b)  and  (d )  have  frequencies  much  above  the 
cutof  f . 

Variations  of  '»  as  a  function  of  path  length  are 
shown  for  rays  it  1  through  it 6.  The  total  path  length 
for  each  ray  is  100  km.  Note  that  rays  ill  and  tt2 
with  f  near  fg=Q  suffer  more  refraction  than  rays 
t)3  through  ttb. 


Figure  4 


Figure  5 


Figure  6 


Figure  7 


Figure  8 


Figure  9 


Variations  of  magnetic  latitude  and  longitude  of 
rays  it 1  through  it 6  are  shown.  The  total  path  length 
for  each  ray  is  400  km.  Air  rays  are  started  at  the 
magnetic  latitude  of  70°  and  longitude  of  0°. 


Growth  and  damping  rates  of  upgoing  AKR  as  a  function 
of  frequency  for  four  wave  normal  angles  are  shown, 
ui  was  computed  by  using  the  exact  distribution 
function  measured  by  the  S3-3  satellite,  run  01. 


The  electron  distribution  function  measured  by  the 
S3- 3  satellite  and  the  loss  cone  boundary  are  shown. 
This  distribution  was  used  in  run  01.  In  runs  02  and 
it  3,  the  region  below  the  loss  cone  boundary  was 
assumed  to  be  empty. 


The  number  of  e-foldings  that  the  wave  electric  field 
grows  or  damps  by  (g),  as  a  function  of  path  length 
(L)  for  run  01. 


Growth  rates  of  AKR  with  M2  at  3  km  below  the 
mirroring  point  Mj  (run  02). 


The  number  of  e-foldings  that  the  wave  electric  field 
grows  by  (g),  as  a  function  of  path  length  (L)  for 
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Figure  10  Growth  rates  of  AKR  with  M2  at  70  km  below  Mj  1  run 

#3)  for  d  between  80°  and  60°. 

Figure  11  Growth  rates  of  AKR  witli  M2  at  70  km  below  Mj  (run 

03)  for  j  between  55°  and  35°. 

Figure  12  The  number  of  e-foldings  that  the  wave  electric  field 

grows  or  damps  by  (g),  as  a  function  of  path  length 
(L)  for  run  #3. 


Figure  13 


Figure  14 


Electron  mean  free  path  plotted  as  a  function  of 
altitude.  In  8a  the  "maximal"  model  for  the 
atmospheric  densities  are  uses,  whereas  in  8b  the 


"minimal"  model  is  utilized. 


Plots  of  r  vs.  pitch  angle  for  electron  energies  of 
0.5  and  5.  keV. 
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